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Second-order post-Hartree-Fock perturbation theory for the electron current 
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Based on the super-fermion representation of quantum kinetic equations we develop nonequilibrium, post- 
Hartree-Fock many-body perturbation theory for the current through a region of interacting electrons. We 
apply the theory to out of equilibrium Anderson model and discuss practical implementation of the approach. 
Our calculations show that nonequilibrium electronic correlations may produce significant quantitative and 
qualitative corrections to mean-field electronic transport properties. We find that the nonequilibrium leads 
to enhancement of electronic correlations. 
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I. INTRODUCTION 

The second order M0ller-Plesset perturbation theory 
(MP2) is a very powerful tool for practical quantum 
chemical calculation.^ MP2 accounts in many cases for 
the large fraction of post-Hartree-Fock electronic correla- 
tions. It gives quite accurate results for intermolecular in- 
teractions, molecular geometries, dissociation energies»i 
In this paper, we propose the extension of the MP2 the- 
ory to nonequilibrium, namely, we develop the second 
order post-Hartree-Fock perturbation theory for the elec- 
tron current through a region of out of equilibrium, in- 
teracting electrons. As an example, we consider electron 
transport through Anderson model, but all our deriva- 
tions are also directly applicable and transferable to a 
general molecular system connected to two metal elec- 
trodes and described by full many-electron Hamiltonian. 

Recently, there have been several attempts to include 
post-Hartree-Fock electronic correlation into transport 
calculations through molecules.^"'' These approaches are 
mostly based on various implementation of nonequilib- 
rium Green's function based GW theory. Our method 
is quite different. We start with kinetic equations for 
many-body density matrix. We convert the kinetic equa- 
tion to the super-Fock space, separate Hartree-Fock and 
the normal ordered parts of the many-particle Liouvil- 
lian, and then develop nonequilibrium many body per- 
turbation theory much along the lines of the traditional 
MP2 theory. 

The rest of the paper is organized as follows. In Sec- 
tion II, we describe quantum kinetic equation and super- 
fermion formalism. Section HI presents the main equa- 
tions of nonequilibrium second-order many-body pertur- 
bation theory. Section IV describes numerical calcula- 
tions. Conclusions are given in Section V. We use natu- 
ral units throughout the paper: fi = ks = |e| = 1, where 
— |e| is the electron charge. 
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II. LINDBLAD KINETIC EQUATION AND 
SUPER-FOCK SPACE 

To be specific in our discussion, let us consider one 
spin-degenerate level (impurity) attached to two macro- 
scopic leads. We begin with the Lindblad kinetic equa- 
tion 



dp{t) 
dt 



aka fi— 1,2 



Pit)}). 



(1) 



Here H is the finite Hamiltonian for the embedded impu- 
rity (the "embedded" here means that the systems con- 
sists of the impurity itself and finite parts of left/right 
leads). The Hamiltonian for the embedded impurity can 
be written in the following form: 

H = Hs + Hb + Ht, 

a 

HT = -tY,Ukaa- + h-c.), (2) 



where , a are fermionic creation and annihilation op- 
erators: a J. (ocr) creates (destroys) an electron with spin 
a in the impurity and a^j,^ (oo-fea) creates (destroys) an 
electron of energy Sk and spin a =t, i in the left {a = L) 
or right (a = R) lead; ria is the number operator for 
electrons of spin a in the impurity; U is the Coulomb 
interaction between two electrons of opposite spin in the 
impurity; t is the tunneling matrix element between the 
states in the leads and in the impurity. For simplicity 
we assume that the tunneling matrix element is real and 
independent of a, k and a. We also assume that the 
left and right leads are identical. Here ^ runs over 2N 

ka 

discrete single particle levels ka coupled directly to the 
Lindblad dissipators. 
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The dissipators have the following form: 

Lakal = y^^kaldaka, Lcka2 = \/'^ka2a\ka^ (3) 
with Tkal = 7fca(l - /fca), ka2 = Ikafka- Here 

fka = [1 + exp[(efea - ^J■a)/Ta]~^ and jka is given by 
the imaginary part of the electrode self energy jka = 
— ■^lm[T,ka]^ The Lamb shift is included into the single- 
particle energt of the leads: Eka — Efc*^^ + Re[SfeQ], where 
^ia bare energy of single-particle levels in the leads. 

We would like to emphasize here that for the steady 
state electron transport calculations the Lindblad master 
equation can be made as accurate and exact as necessary 
by the increasing the number of lead states k included 
into the Hamiltonian for the embedded impurity. So the 
master equation ([T|) can be considered as the exact start- 
ing point for first principles electronic transport calcula- 
tions. 

Let us begin to work with the kinetic equation ([T|). 
First, we re-write it in the super-fermion representation. 
For many-particle quantum systems the density matrix 
p(t) and the Hamiltonian are operators in the Fock space. 
The Fock space can be defined by some orthonormal com- 
plete set of basis vectors: 

^|n)(n|=/, {n\m) = Snm- (4) 

n 

Let us introduce the additional Fock space which is iden- 
tical copy of the initial Fock space 

^\n){n\^I, {n\m) ^ 5nm- (5) 

n 

We denote all vectors and operators in this additional 
Fock space by "tilde". The vectors \n) and \n) span the 
so-called super-Fock space, which is a direct product of 
the original and the "tilde" Fock spaces. In the super- 
Fock space we define "left vacuum vector" 

|/) ^^a„|n)«)|H}, (6) 

n 

and " nonequilibrium wavefunction" : 

\p{t))=p{t)\I), (7) 

where a„ (lanl = 1) is an arbitrary phase factor. Now 
the nonequilibrium average can be written as 

{A{t))=i:T[Ap{t)] = {I\A\p{t)). (8) 

We consider a system which consists of fermions dis- 
tributed over N levels. If we take vector \n) and \n) to 
be the particle number eigenstate, |n) = \n1n2 ■ ■ - nN) , 
n) = |nin2 . . . riAr), and choose q;„ = (— i)"i+"2+---+niv^ 
then we can readily demonstrate by the straightforward 
algebraic manipulations that 

a,\I) = -ia]\I), a]\I) = -ili,\I). (9) 



This is so-called "tilde conjugation rule". It transforms 
the original operators to the tilde operators and it is one 
of the most important relations. In particular, it follows 
from ^ that the left vacuum vector, (/| , is always the 
vacuum for a}- — iuj and -I- iaj operators. Moreover, it 
shows that, in some sense, the particle creation (annihila- 
tion) in ordinary Fock space is equivalent to the particle 
annihilation (creation) in the "tilde" space. 

We can formulate the following rules for the operations 
in the super-Fock spaces^"— 

1. The left vacuum vector |/) and the nonequilibrium 
wavefunction \p(t)) = p{t)\I) are invariant under 

the tilde conjugation \p{t)) — \p{t)), \I) — \I), and 
{I\p{t)) = l. 

2. Tilde conjugation rules: The main rule is aj\I) = 
—ia^I), — —iall) and as consequence the 
double tilde conjugation does not change the oper- 
ator A and {aA + C2B) = c^A + c^B. 

3. Evolution of the system is described by the 
time-dependent Schrodinger equation i-^\p{t)) — 
L|p(t)), where the Liouvillian is obtained from the 
corresponding master equation for the density ma- 
trix by the tilde conjugation rules. The nonequi- 
librium average is given by {A{t)) — {I\A\p{t)) and 
{I\L = 0. 

If we act by the Lindblad equation on the "left vac- 
uum vector" |/), employ the tilde conjugation rules and 
use the fact that the density matrix p = p{a\a) is the op- 
erator in original Fock space therefore it commutes with 
all "tilde" operators, then the Lindblad master equation 
becomes the time-dependent Schrodinger equation in the 
super-Fock space 

^±\pit)) ^ L\p{t)) (10) 

where the Liouvillian L is given by 

L = Ls + Lb + Lt, (11) 

where = Hs — Hs and Lt = Ht — Ht are both 
Hermitian, while 

Lb = Hb - Hb - lY^n^ka (12) 

aka 

includes the non-Hermitian part which is responsible for 
the dissipation in the system: 

— kaiaakaaaka + ^ ka2a\ka'^'^aka) + '^^ka2- 

(13) 

In the above equations all tilde operators are obtained 
from non-tilde ones by the tilde conjugation rules. 
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In this paper, we focus on nonequilibrium steady state, 
where the density matrix \p{t)) has aheady reached its 
asymptotic steady state |poo) and does not anymore de- 
pend on time. Therefore, the electron transport prob- 
lem is reduced to the problem of finding the eigenvector 
with zero eigenvalue of complex, non-Hermitian, finite- 
dimensional Liouville operator (jlip acting in the super- 
Fock space 

L\poo)=0- (14) 

Once the steady state density matrix is found, we can 
compute the current as 

Ja = -it^(/| {alka^a - ala^ka) IPoo)- (15) 

trk 

In the next section, we show how one can find \poo) per- 
turbatively beyond nonequilibrium Hartree-Fock approx- 
imation. 



III. NONEQUILIBRIUM MANY-BODY 
PERTURBATION THEORY 

Before we start to develop many-body perturbation 
theory, we make the important remark on the nota- 
tion abuse in this section of the paper: only cre- 
ation/annihilation operators written with letter a (such 
as for example aka and ci\^) are related to each other by 
the hermitian conjugation; all other creation and anni- 
hilation c operators are "canonically conjugated" to each 
other, i.e., for example, does not mean (c)^ although 
{c,ct}^l. 

Let us first perform the normal ordering of the Liou- 
villian, and find the nonequilibrium Hartree-Fock density 
matrix \p^) and corresponding nonequilibrium quasi- 
particles (i.e. quasiparticles, which have the Hartree- 
Fock density matrix as a vacuum in the super-Fock 
space). 

Using the Wick theorem we rewrite the Liouvillian pT|) 

as 

i: = L(")+L', (16) 

where L' is the normal ordered part of the Liouvillian 

L' = U : {a\a-^a\ai - t.c.) :, (17) 

and the notation (t.c.) means the tilde conjugation (i.e. 
Oct — ^ Co-, etc.). The normal ordering is asymmetric: it 
is performed with respect to (/| from the left and the 
nonequilibrium Hartree-Fock density matrix from 
the right. 

The mean-field part of L is 

L^"") = Lf + Lt + Lb, (18) 



where Lg is a quadratic part of Ls in the Hartree-Fock 
approximation 

= E(^0 + C/(n-a))(4a. - a^fla). (19) 
a 

The mean-field population {n„) = (/| no- Ipm'*) is spin in- 
dependent and {a\a-a) ~ due to the symmetry of the 
problem. The quadratic part of the Liouvillian, L*^°\ 
can be diagonalized exactly in terms of nonequilibrium 
quasiparticle creation and annihilation operatorsJ^ As a 
result we get 

L(") = ^(f^„cLc.„ - f^:4„c.„), (20) 

na 

where the nonequilibrium quasiparticles are defined as 

kOL 

Can = V'nacr + i^Pn{a\ + la^) 

+ '^{'4^n,kaaaka + i^n,ka(a\ka + io,ka)) , (21) 
ka 

and cj.„ , Can are obtained from , Can by the tilde con- 
jugation rules. The operators c^, c, 7?, c are connected 
to a\ a, a\ a by canonical (but not unitary) transfor- 
mations. Nonequilibrium quasiparticle creation and an- 
nihilation operators obey the fermionic anticommutation 
relations. By the construction, (/| is the left vacuum for 
cjj, cjj operators. The vacuum state for c„, c„ operators, 
\poo), is automatically the zero-eigenvalue eigenstate of 
the mean-field Liouvillian L^^\ and therefore, it is the 
steady state density matrix in the Hartree-Fock approx- 
imation: 

L(°)|p£))=0, (/|p(S)) = L (22) 

In other words, the mean-field density matrix is the den- 
sity matrix which does not contain nonequilibrium quasi- 
particle excitations. 

The amplitudes ipn, 4'n,ka and quasiparticle energies 
^n {n = 1,...2N + 1) are the solution of the follow- 
ing eigenvalue problem (obtained from the equations-of- 
motion [L, 4„] = 1^„4„) 

ka 

Ekai^n.ka — t'ipn = ^ni^n.ka, (23) 

where Eka — £ka — ilka and Ehf = + U{na) is the 
Hartree-Fock single particle energy. The amplitudes ipn, 
ipn.ka are normalized as 

Ipn-lpn' +^^i>n,kai>n',ka = 5nn'- (24) 
ka 

The amplitudes ipm fn.ka satisfy the following nonho- 
mogeneous system of linear equations (obtained from 
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ka 



ka 



i^ka ~ ^n)fn,ka — t^n — —tfkai^n- 



(25) 



This is a set of nonlinear equations since the Hartree-Fock 
single particle energy shf depends on amplitudes 



(26) 



which in turn depends on Ehf- We solve these equa- 
tions numerically by self-consistent iterations. We first 
guess the population (ria-), then for a given ehf we solve 
the eigenvalue problem (1^51) . Then we solve the linear 
system of equations (1251) with the known quasiparticle 
spectrum i7„ and known amplitudes V'n, ^n,ka- Then we 
compute new nonequilibrium mean-field population {n„) 
of the impurity and continue this loop for self-consistent 
iteration until the variation in the impurity population 
between the iterations becomes negligible. 

Now we are in a position to develop the post-Hartree- 
Fock many-body perturbation theory, which takes into 
account the normal ordered part of the Liouvillian (IT71) . 
We introduce the continuous real parameter A, which will 
be set to unity in the end of the calculations 



L = + XL'. 



(27) 



The normal ordered post-Hartree-Fock perturbation L' 
in the nonequilibrium quasiparticle basis is given in the 
Appendix. We expand the exact steady state density 
matrix in powers of A, 



Poo)= |p^^) + A|p^^)+A^|p^O + 
Sp)\ 



(28) 



where (/| = 0, p > 1. If we require that \poo) is the 
steady state density matrix for the full Liouvillian (ITil) . 
we get the following recurrent relations for the pth-order 
correction to the Hartree-Fock density matrix: 



(29) 



In the case of the second-order perturbation theory this 
recurrent relation simply becomes the systems of two 
equations 



(30) 
(31) 



Before solving Eas. ipOl I5T]) for and \p^^) let us 

represent the current p5p in terms of nonequilibrium 
quasiparticles and understand which configurations in 
nonequilibrium density matrix give nonzero contribution 
to the expectation value of the current. Using the trans- 
formation inverse to (1211). we find 



Ja=-it 22 {(C,fcaVm " ^nV*m,ka){I\ Cna^^a \Poo) 
aknni 

IPoo)}- (32) 



The first term here is fully accounted within the Hartree- 
Fock approximation 



kn 



(33) 



and does not contribute to the expectation value over 
the correlated density matrix. The pth-order perturba- 
tive theory correction to the HF current comes from the 
second part of Eq. 



fcmn 

= ~mmJ2rn,kc.^mF^l (34) 



kran 



where FmX is the expansion coefficient in 

\FoD I ' / , ^ mn^am^an Iroo / ^ 



(35) 



and F^t = {F^IY as follows from \pt^) - 

Thus, to find the pth-order correction to the Hartree- 
Fock current we need to calculate Fmn- Since L' is normal 
ordered it does not contain terms quadratic in nonequi- 
librium quasiparticle creation operators. Therefore, as it 
can be easily seen from Eq. ipO)) . Fmn — 0, and the first 

(2) 

nonvanishing correction to the current is JA . 

Eq. ([5n)l gives the first-order perturbation theory cor- 
rection to the Hartree-Fock steady state density matrix. 
It contains the mixture of four nonequilibrium quasipar- 
ticle excitations: 



(1) 



klmn 



(36) 



where 



G 



(1) ^klmn 



(37) 
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and G^;)„„ = G^^^^^,. Substituting ^ into (EH) we ob- 



tain 



where 



|P£^> = * E P^nlclj^an \pt^) + ■ ■ • , (38) 



p(2) ^ 1 \^ f r^'^) r*^^ _ 7-*(5) ^*(i) -j. 

^ mn Q _ O* / ' \ milk kmi nijk kjmil' y"^^ ) 



ijk 



Inserting (p9)) into (p4)) we get the second-order perturba- 
tion theory correction to the Hartree-Fock current ([55]) . 

The omitted terms in ([551) contain four, six, and eight 
excited nonequilibrium quasiparticles and these configu- 
rations do not contribute to the current in the second- 
order perturbation theory. However, they contribute in 
higher orders. We also note that each term in \p'£^) 
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contains even number of excited quasiparticles and the 
number of tilde and non-tilde creation operators in the 
excited quasiparticle configurations is the same in all or- 
ders of perturbation theory. 

To this point, the calculations are valid for arbitrary 
left and right leads. For the identical leads, Eqs. (|33l34p 
can be further simplified. Indeed, if we have identical left 
and right leads, i.e., if EkL = EkB. = Ek, Eq. (|23|) has 
N eigenvalues given by iln = Ei. For these eigenvalues 
i>„ = 0, and ^|:n,kL = -i'n.kR = -4z2^^/2. The re- 
maining A^-l- 1 eigenstates have 4'n,kL = 4'n,kR and hence 
they do not contribute to the sum over n in Eqs. (|33|34|) . 
Thus, we obtain 



(40) 
(41) 



where the upper (lower) sign corresponds to a — L (R) 

r{5) 



("5) 

In addition, because of ■0n = 0, i„/,i, — and the second 



term in Eq. ((39l) vanishes. 

Moreover, using explicit expression for tpn in the case 
of 0„ = El, 



2-1/2. 



tiflL - flR) 



(£HF-i?/) + 2i2^-^— - 
k [El - E^ ) 



(42) 



we can readily demonstrate that if the macroscopically 
large part of the leads (A^ ^ oo, 7^ — >• ) is included 
into the embedded Hamiltonian ([2]) , Eq. (PO]) becomes the 



Landauer formula for the current in the Hartree-Fock 
approximation (seei^ for details). 



IV. NUMERICAL RESULTS 

We perform numerical calculations assuming that each 
lead has N discrete, spin degenerate energy levels uni- 
formly distributed in the bandwidth [Emm ■ Emax] — 
[—5 : 5]. The tunneling coupling strength t is computed 



from the F — 2iTrit^ = 1, where 77 N/{E„ 



Emin ) 



is the density of states in the leads. In all our calcu- 
lations we take N = 1600 and 7 = 2A£, where Ae is 
the energy spacing between states in the leads. Under 
this choice of the parameters the imaginary part of the 
self-energy becomes energy independent constant F = 1 
within the bandwidth of the leads and vanishes outside 
the bandwidth. We use a symmetric applied voltage: 

^JiLM = ±0.5^. 

To verify that the sufficient number of the leads levels 
is included into the embedded Hamiltonian ([2]) we calcu- 
lated the Hartree-Fock current for two different regimes, 
for fixed applied voltage V and varying central cite en- 
ergy Eq and vice versa, and then compared our results 
with the Landauer, i.e. exact, mean-field current. We 
found the disagreement between obtained results is neg- 
ligible (the maximum deviation < 1% of the value of the 
current). 

Fig. [1] shows the current plotted as a function of the 
impurity level for different values of the Coulomb inter- 
action energy U . As we can see the HE current jf^ 
becomes broader with increasing of U and reaches its 
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maximum value at the symmetric point eo = —U/2. The 
second order perturbation theory correction to the cur- 
rent is symmetric with respect to — —U /2 and can be 
both positive and negative. We have numericahy found 

(2) 

that Jj^ < only when the Hartree-Fock level Ehf is 
between and ^r. For C/ < 1 the second order pertur- 
bation theory correction to the current is very small, that 
means the Hartree-Fock approximation accounts for the 
main part of the electronic interactions. For larger U the 
second order perturbative corrections change the behav- 
ior of the current both qualitatively and quantitatively. 
For example, when [/ = 3 and V = \ the electronic cor- 
relations produce the dip in the current in the vicinity of 
the symmetric point eo = —U/2. 

We found that the larger the voltage (i.e. the further 
away we are from the equilibrium), the more important 
role the nonequilibrium electronic correlations play. This 
is evident from the comparison of the left (computed at 
V — 1) and right (computed a,iV = 0.1) panels of Fig. [T] 
It means that the nonequilibrium dynamics and elec- 
tronic correlations are entangled in nontrivial way and 
the nonequilibrium enhances and amplifies the role of 
electronic correlations. 



V. CONCLUSIONS 

In this paper, we developed the second-order post- 
Hartree-Fock perturbation theory for the electron current 
through a region of out of equilibrium, interacting elec- 
trons. As an example, we considered electron transport 
through out of equilibrium Anderson model, although 
all our derivations are also directly applied to a gen- 
eral molecular junctions described by full many-electron 
Hamiltonian. We started with the Lindblad kinetic equa- 
tion for the embedded molecular system. We converted 
the kinetic equation to super-fermion representation and 
define nonequilibrium quasiparticles within Hartree-Fock 
approximation. Then we applied the Wick theorem and 
perform the normal ordering of the Liouvillian with re- 
spect to the vacuum for nonequilibrium quasiparticles. 
We developed the second-order post-Hartree-Fock per- 
turbation theory by admixing two- and four-quasiparticle 
configurations to the nonequilibrium vacuum. Our nu- 
merical calculations for out of equilibrium Anderson im- 
purity demonstrated that nonequilibrium electronic cor- 
relations may produce significant quantitative and qual- 
itative corrections to Hartree-Fock electronic transport 
properties. We also found that the nonequilibrium en- 
hances the role of electronic correlations. 
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Appendix A: Normal ordered Liouvillian in the basis of 
nonequilibrium Hartree-Fock quasiparticles 

Here we give the explicit expression of L' in terms 
of nonequilibrium quasiparticle creation and annihilation 
operators. Using the inverse transformation 

71 n 

aa ='^i^nCan~i'^'PnCln (^-1) 
n n 

we obtain 

L' —U : (a|a-|-a|aj_ — t.c.) : 

— X {(-^fe/rnn'^fe-fC|^Cmj^C„^ — t.C.) 
klmn 

, r (2) J J p p 
+ ^''klmni'^Xf^l^'^miCn^ + c\^c\^Cm^Crn_ 

(A2) 

where l[,',^„„ are given by 

^kllnn =U{'lpk1pl - ll>kVl - VfeV'OV'mV'n, 
^kLn =U[{lpk - Vk){'4^l - ^l)^*m^n 

^klmn =U{ipkll'! - TpkViW^i^n, 
LkLn = - U[{ijk - Vk){i^l ~ (pi)V*m 
+ Vkfli^rn - 'P*m)]i^n, 

It is notable that because oi {I\ L ~ 0, L' does not contain 
terms of four annihilation operators. 
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